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We calculate the energy spectrum of closed strings in SU (N) gauge theories with N = 2,3,4,6,8 
in 2 + 1 dimensions to a high accuracy. We attempt to control all systematic errors, and this allows 
us to perform a precise comparison with different theoretical predictions. 

When we study the dependence of the string mass on its length L we find that the Nambu-Goto 
prediction is a very good approximation down to relatively short lengths, where the Liischer term 
alone is insufficient. We then isolate the corrections to the Liischer term, and compare them to 
recent theoretical predictions, which indeed seem to be mildly preferred by the data. 
When we take these corrections into account and extract string tensions from the string masses, 
we find that their continuum limit is lower by 2% — 1 % from the predictions of Karabli, Kim, and 
Nair. The discrepancy decreases with N, but when we extrapolate our results to N = °° we still 
find a discrepancy of 0.88% which is a 4.5(7 effect. 
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1. Introduction 

There are many indications for a connection between QCD and string theory, which appear in 
a wide range of experimental phenomena and theoretical works [fflj. Here we are interested on the 
QCD side and study the energy spectrum of the 'original' string - the flux tube of pure Yang-Mills - 
with an emphasis on its large-A^ limit. We focus on three Euclidean dimensions, and our motivation 
is two-fold. First, we wish to perform a precise test of the remarkable work of Karabali, Kim, and 
Nair [§]. This work yields the following prediction for the string tension a 



^v-V^^' (L1) 

where g is the super-renormalised coupling 1 . A previous comparison of Eq. ( \l.l\ > to lattice results 
for N = 2, 3,4, 6 [§| showed that Eq. ( [01 ) is higher from the data by about 2% - 1%, but also that 
this discrepancy decreases with N. This and the fact that [§] predicts no screening of zero Af-ality 
charges suggests the possibility that the analysis in [^] may be exact at N = °°. This would appear to 
be contradicted by the N — > oo extrapolation of but the results presented there included several 
unestimated systematic errors which, while small, may be significant at the 1% level. Two of these 
lead to an underestimate of the lattice string tensions, and we check whether by controlling them, 
the lattice results become consistent with Eq. ( |1.1[ ) at = oo. 

The first error rises in the process of extracting the string tension from the string mass, where 
it is typical to neglect corrections that are sub-leading to the Liischer term. Indeed, as the string 
length L becomes larger, these corrections decrease as l/L , but we wish to estimate them on a 
quantitative level. The way the string mass changes with L is also interesting theoretically (and 
this constitutes our second motivation for this study). It gives us information on the effective string 
theory of the flux tube, that can be compared with older [Q] and more recent [j|] predictions. The 
second systematic error that we remove is the neglect of 0(g 4 ) terms in the continuum extrapolation 
of the string tensions. 

Here we present an analysis of these systematic errors, and in the case of the 0(1/L 4 ) error, 
we compare to the predictions of After removing both systematics we compare the resulting 
string tensions to Eq. The results presented here were obtained in a preliminary analysis 

of the data. Publications with a more extensive analysis, that will also include the raw data, are 
forthcoming (Bp. 



2. Methodology 

There exists a plethora of lattice works studying properties of flux-tubes by looking at open 
or closed strings (for example see the review in [§]). To avoid 'contamination' from perturbative 
effects, and to compare directly with the spectrum of the effective string theories, we choose to 
measure the large distance exponential fall-off of correlation functions of strings closed around the 
spatial torus. Our analysis proceeds in two stages. 

'Recall that in 2 + 1 dimensions the coupling g 2 has dimensions of mass. 
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2.1 Stage I - the string mass dependence on the string length 

We first work with a fixed lattice spacing (by fixing the lattice coupling) and study the depen- 
dence of the string mass m on its length L. Fitting these results we are able to test the theoretical 
predictions in [Q, gj. In addition, these fits provide us with a practical tool to extract string tensions 
from strings which are shorter then usual, and with a quantitative estimate of the systematic error 
induced by using the Liischer term alone. 

2.2 Stage II - continuum extrapolation of string tensions 

Here we perform measurements of string masses with a fixed physical length L and different 
lattice spacings. We choose L\fc > 3 — 3.5 and use the fits we obtain from stage I to extract string 



tensions. We then extrapolate these to the continuum limit and compare with Eq. ( |1.1| ) for all values 
of N as well as for the extrapolation to N = °°. 

2.3 Lattice construction 

We define the gauge theory on a discretized periodic Euclidean three dimensional space-time 
with A^o x Ni x N2 sites, and perform Monte-Carlo simulations of a simple Wilson action. We use 
the Kennedy-Pendelton heat bath algorithm for the link updates, followed by five over-relaxations 
of all the SU(2) subgroups of SU(N). We measure correlation functions of Polyskov loops that 
wind around the 6 direction so L = ciNq. The correlations are measured along direction 1 and we 
project to zero transverse momentum by averaging over direction 2. Although here we mainly 
present results for the ground state, we study all possible values of the Af-ality k = 1,2, . . . , [N /2], 
and their first few excited states. Using improved operators []7|] we are able to obtain overlaps 
which are almost perfect for the ground state of k = 1, but somewhat lower for the excited states. 
To avoid finite volume effects we increase the lattice in the 1 and 2 directions for the shorter strings 
[@]. Stage-I is studied for N = 3,4,6,8, and 1.3-1.6 < Ly/a < 3-6.2, while in stage-II we study 
SU (2) as well and restrict to La/ct > 3 — 3.5 with 0.1 < a^fo < 0.75. The raw data will be presented 
in 



3. Results - dependence of the string mass on its length 

We first check the universality class of the string by looking at the behaviour of an effective 
central charge defined as = a — KC f^ ■ Here C e e (L) should be 1 at L 3> 1/ \fc for the bosonic 
string, and is obtained from our data by performing fits to successive pairs of adjacent points. The 
results are shown in Fig. [j] where we see that for yfoL > 3 the charge C e ff becomes consistent with 
1. 

Next we fit our data to the general form 2 

^ = rT _A__^ D E 

L L 2 (oL 2 )L 2 (oL 2 )V 2 L 2 (oL 2 ) 2 L 2 ' K ' 

Here we focus on the following fits. First we fix the 0(1 /L 2 ) to be the Liischer term with B = n/6 
and (i) let C be a free fit parameter but fix D = E = (ii) follow the theoretical predictions fj3j] and 

2 The 1/L 3 term is missing since it was unacceptable fit for our data, and was suggested by Liischer and Wiesz to be 
disfavoured theoretically o]. 
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Figure 1: Effective central charge as a function of the string length, for SU (3),j8 = 14.7172 (blue circles), 
SU(4),P =28.00 (red dots), SU(6),P = 59.40,90.00 (blackstars and squares), andSI/(8),j3 = 108.00 (pink 
crosses). 



constrain C = 7T 2 /72 while fitting with D ^ 0,E = or D = 0,E ^ 0. We also compare our fits to 
what one obtains when one uses the Liischer term alone m Luscher (£) /L = G — -Mj, or the Nambu- 
Goto formula m^Q{L)/L = O^Jl — j^jj as obtained by Arvis [Q]. (Here a is the string tension 
obtained from our fit) As a demonstration, we present results for SU (6) and /3 = 90.00 in the left 
panel of Fig. ||, where one can see that the Liischer term is a good approximation for yfoL > 3 while 
the Nambu-Goto prediction works remarkably well down to almost the deconfinement length, but 
not exactly. 
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Figure 2: Left: m(L)/L for the ground state closed string of SU (6) at j3 = 90.00. The fit (red line) gives 
ci\fa = 0.172148(70). Substituting this in rauischei- an d m NG gives the magenta and blue lines. Right: m(L)/L 
for the first excited state of the k = 1 sning. The lines are plots of mmscher an d m NG for a quantum n = 1 and 
<7 as extracted from the ground state. 

The results for all the other gauge groups are presented in Table [l], where one can see that the 
scaling and l/N effects are relatively small. Our best fits in the case (i) and (ii) have comparable 
confidence levels of 43% — 90%, but a mild preference is seen toward option (ii). In apparent 
contrast to that, the coefficient C obtained in the fits of type (i) turns out to be about twice as 
large as predicted in Nonetheless, when we perform fits in which both C and either E are free 
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(denoted as type (iii) in the last two columns in Table [T]) we find that the data is almost consistent 
with C ~ 7T 2 /72 as predicts. 



Data set 


fit (i), D = E = 


fit (ii), C = K 2 /12 


fit (iii), D = 


N 


P 




D E 


c p 


3 


14.7172 


0.261157(82) 1.96(8) 


0.230(6) 0.299(8) 


1.32(12) 0.23(3) 


4 


28.00 


0.25198(16) 2.26(8) 


0.23(4) 0.36(2) 


1.36(32) 0.25(9) 


6 
6 


59.40 
90.00 


0.27870(12) 1.91(6) 
0.172148(70) 1.85(6) 


0.17(1) 0.25(2) 
0.16(1) 0.22(2) 


1.45(23) 0.12(6) 
1.21(24) 0.16(6) 


8 


108.00 


0.27402(34) 1.58(2) 


0.13(4) 0.19(6) 


0.77(32) 0.28(9) 



Table 1: The results of our fits of type i-iii (see text). The values of ay/a for fits of type ii,iii were consistent 
mostly within 1 .5(7 with those for type i, which are presented in the table. 

Finally we can use the value of a obtained from the ground state of the string, and substitute 
it in the Liischer and Nambu-Goto formulas, but with a quantum n = 1 (see for example the work 
of Arvis [Q]). We compare the result to the mass of the first excited state that we measure. The 
comparison is presented in the right panel of Fig. |2[ We see that the Liischer formula does not 
describe the data at all, but that, surprisingly, the full Nambu-Goto describes the data quite well. 
This seems to suggest that it is more natural to fit the data for m 2 rather than for m, since the former 
has the full Nambu-Goto expression as a zeroth approximation. An analysis in this form will be 
presented in a forthcoming publication [Sj. 

To conclude this section we find that the systematic error induced by assuming C = D = E = 
for strings with \J~cL ~ 3 is to underestimate the string mass by about a third of a percent. In the 
next section we use this result and extract string tensions using the form Eq. (3.1) with B = n/6, 
D = E = and C as given by the fits to the different gauge groups (we assume that the scaling 
violations in C are small, as observed for SU (6) - see Table fil). 



4. Results - extrapolation of string tensions to the continuum 

Equipped with the fitting formula for m(L) we extract string tensions from a series of mass 



measurements where we keep the physical length of the string to be Lyff > 3 — 3.5 and change the 
lattice spacing by changing j8 . We perform a continuum extrapolation with the ansatz 

bo b\ 



->MF 



2N 2 



g 2 N 



■>MF 



A 2 



(4.1) 



MF 



Here, Pmf is the mean field improved coupling Pmf = P • («p). for which we measure the expecta- 
tion value of the lattice plaquettes (u p ) at each /3. Fitting with Eq. ( JO| ) is the way we remove the 
systematic error related to the 0(g 4 ) terms that we mention in the introduction. As a demonstration 
we present the results for SU (4) and SU (6) in the left panel of Fig. || where it can seen that b\ > 
and therefore that the removal of this systematic error increases the estimate of the continuum 
string tension. Nonetheless, as seen in the Figure, the extrapolation of the string tension are still 
lower compared to the value predicted by Karabali, Kim, and Nair (KKN) for these gauge groups. 
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Figure 3: Left: The dimensionless quantity )3mf^t as a function of the improved coupling I/Pmf for 
St/ (4) (lower plot, in red) and St/ (6) (higher plot, in blue). The error bars at 1 /Pmf = denote the result 
of the continuum extrapolation, while the horizontal bars denote the values predicted by Karabali, Kim, 
and Nair [||]. Right: The ratio between the prediction of [Q] and our data, r, as a function of l/N 2 , for 
iV = 3,4,6,8. The error bar at 1 /N 2 = denotes the extrapolation to the N = °° limit, while the horizontal 
line at r = 1 denotes where the prediction of [||] exactly matches the lattice result. 

In the right panel of Fig. ^] we present the ratio r = £7=71 ^kkn between the Karabai-Kim-Nair 

(V IS J Lattice ._, 

prediction [g] and the lattice data (the numerical values are given in Table |2|). 

5. Summary 

We measure string masses in SU (N) gauge theories with Af = 2,3,4,6,8 in 2 + 1 dimensions 
using lattice techniques. Here we present the results of two studies performed with these masses. 
In the first we investigate the behaviour of the closed string masses as a function of their length at a 
fixed lattice spacing. We find that for ^foL > 3 our data is consistent with a Liischer term of a unit 
central charge. For shorter strings the Liischer term is insufficient and we find that the Nambu-Goto 
prediction works surprisingly well, but is not exact. We fit the deviations from the Liischer term 
and compare the result to the recent theoretical predictions which seem to be mildly supported 
by the data. We estimate that the systematic error induced by neglecting these sub-leading terms at 
L\fo ~ 3 is to underestimate the string masses by about a third of a percent. 

In the second study we calculate string tensions by restricting to string lengths of \fcL > 
3 — 3.5 and change the lattice spacing. Using the fit from the first study we extract the string 
tensions and extrapolate to the continuum limit. The results are given in Table |2[ where we also 
give the Karabali-Kim-Nair (KKN) predictions. As indicated in the table, when we extrapolate to 
N = 00 we find that the lattice result is lower by about 0.88% then the prediction of [Q] which is a 
4.5a effect. 

6. Future prospects 

An additional systematic error is related to the contamination from excited states, which we 
take partially into account by fitting correlation functions with a single exponential at large t. Elim- 
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N = 2 


7Y = 3 


N = 4 


N = 6 


N = 8 


N = oo 


KKN 


0.17275 


0.1881 


0.1931 


0.1967 


0.19791 


0.19947 


Lattice 


0.16678(42) 


0.18425(28) 


0.1898(4) 


0.1946(4) 


0.19580(47) 


0.1977(4) 



Table 2: k = 1 string tensions for SU(N) pure gauge theories in the continuum and the predictions of 
Karabali, Kim, and Nair (KKN) [@. 



inating this effect will tend to give a lower string tension, which would increase the discrepancy 
with Eq. ([fl]). Also, since in the case of /^-strings our operators have typically a lower overlap than 
in the case of k = 1, it is possible that one will obtain a significantly lower ^-string tension there. 
This may bring the results numerically much closer to Casimir scaling, which may have interesting 
theoretical implications, and an initial analysis of this issue will be present in [fl]. Other directions 
of research include a detailed spectrum calculation of excited string states and of glueballs [^] with 
a larger basis of operators that could be compared with recent theoretical predictions [ |10| ] 
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